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Abstract 

Given a set of integers W, the Partition problem determines whether W can be divided into 
two disjoint subsets with equal sums. We model the Partition problem as a system of polyno- 
mial equations, and then investigate the complexity of a Hilbert's Nullstellensatz refutation, or 
certificate, that a given set of integers is not partitionable. We provide an explicit construction of 
a minimum-degree certificate, and then demonstrate that the Partition problem is equivalent 
to the determinant of a carefully constructed matrix called the partition matrix. In particu- 
lar, we show that the determinant of the partition matrix is a polynomial that factors into an 
iteration over all possible partitions of W. 
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1. Introduction 

The NP-complete problem Partition |9|) is the question of deciding whether or not 
a given set of integers W = {wi, . . . ,w n } can be broken into two sets, / and W \ I, 
such that the sums of the two sets are equal, or that J2wei w = YlweW\i w - Since it is 
widely believed that NP ^ coNP, it is interesting to study various types of refutations, 
or certificates for the non-existence of a partition in a given set W. 

In this paper, we study the certificates provided by Hilbert's Nullstellensatz (see (jU 
[SJ [SJ Uni H2J) and references therein) . Given an algebraically-closed field K and a set of 
polynomials /i, . . . , f s £ K[xi, . . . , x n ], Hilbert's Nullstellensatz states that the system 
of polynomial equations /j = f 2 = ■ ■ ■ = f s = has no solution if and only if there 
exist polynomials j3\, . . . , j3 s £ K[xi, . . . , x n ) such that 1 = Pifi ■ We measure the 

complexity of a given certificate in terms of the size of the (3 coefficients, since these arc 
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the unknowns we must discover in order to demonstrate the non-existence of a solution 
to /i = /2 = • • • = f s = 0. Thus, we measure the degree of a Nullstellensatz certificate 
as d = max{deg(^i), . . . , deg(/3 s )}. 

There is a well-known connection between Hilbert's Nullstellensatz and a particular 
sequence of linear algebra computations. These sequences have been studied from both 
a theoretical perspective (jHIHl), and a computational perspective ([7Jin|). In Q, Buss and 
Pitassi show that a polynomial system loosely based upon the "pigeon-hole principle" 
requires a [log tt,J — 1 Nullstellensatz degree certificate. However, when the system of 
polynomial equations f\ , . . . , f s models an NP-complete problem, the degree d is likely 
to grow at least linearly with the size of the underlying NP-complete instance ([TT|). In 
other words, as long as P ^ NP, the certificates should be hard to find (i.e., the size of 
the linear systems involved should be exponential in the size of the underlying instance) , 
and as long as NP ^ coNP, the certificates should be hard to verify (i.e., the certificates 
should contain an exponential number of monomials). 

For example, consider the NP-complete problem of finding an independent set of size 
A; in a graph G. Recall that an independent set is a set of pairwise non-adjacent vertices. 
This problem was modeled by Lovasz (|10|) as a system of polynomial equations as follows: 

x\ — Xi = , for every vertex i G V(G) , 
XiXj = , for every edge £ E{G) , 

-k + J27=i Xi=Q . 

Clearly, this system of polynomial equations has a solution if and only the underlying 
graph G has an independent of size k. For example, consider the Turan graph T(5, 3): 

3 4 




By inspection, we see that size of the largest independent set in T(5, 3) is two. Therefore, 
there is no independent set of size three, and using the connection between Hilbert's 
Nullstellensatz and linear algebra (described more thoroughly in Sec 3), the authors of 
(JSJ) produce the following certificate: 

1 = f — i (xiX2 + X 3 Xa) - g (x\ + X 2 + X 3 + X 4 + IE 5 ) - 1^ (X 1 + X2 + X S + X4 + x 5 - 3) + 
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(l Xi+ 3 X2+ 3) XlX3+ (l X2 
(J^Jx 2 X4 + (J^jx 2 x 5 + Qx 4 + ~^x 3 x 5 + Qjaysg + Qz2 + gj( a l _ + 

(j$ Xl + e) ^ ~ X2 ^ + {^ X4 + ^ ~ X3 - ) + + e) ^ ~ X ^ + (s)^ x * ~ X ^ ' 

The combinatorial interpretation of this algebraic identity is unexpectedly clear: the size 
of the largest independent set is the degree of the Nullstellensatz certificate (i.e., the 
largest monomial X1X2 corresponds to the maximum independent set formed by vertices 
{1, 2}), and the coefficient fi\ contains one monomial for each independent set in G. The 
combinatorial interpretation of these certificates is proven in ([5| by De Loera et al. only 
in terms of monomials: the relationship between the numbers such as 1/3 and 1/6 and 
the independent sets of the underlying graph is not clear. 
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In this paper, we model the Partition problem as a system of polynomial equations, 
and then present a combinatorial interpretation of an associated minimum-degree Null- 
stcllensatz certificate. However, the focus of our combinatorial interpretation is not only 
on the relationship between partitions and monomials, but also on the relationship be- 
tween partitions and numeric coefficients (i.e., the numbers 1/3 and 1/6). In Section 2, 
we present an algebraic model of the partition problem and describe a minimum-degree 
Nullstcllcnsatz certificate. In Section 3, we describe the connection between Hilbert's 
Nullstellcnsatz and linear algebra, leading to the construction of a square system of lin- 
ear equations, forming what we call the partition matrix. In Section 4, we prove our main 
result: the determinant of the partition matrix represents a brute-force iteration over all 
the possible partitions of the set W, a polynomial we refer to as the partition polynomial. 

We conclude our introduction with an example. Let W — {wi,w 2 ,W3, w 4 } 7 and we see 
that the determinant of associated partition matrix is as follows: 



det 



U 



W4 w 3 w 2 w\ 

(1)3 1«4 U)2 W\ 
W2 W4 W3 W\ 

iti W4 W3 W2 
W2 w 3 W4 t»i 
w\ W3 W4 W2 
wi W2 W4 W3 

1»i 11)2 t»3 !»4 



\ 



(wi + ui 2 + w s + W4)(-Wl + U>2 + M>3 + UI4) 
(wi - u>2 + w s + W4)(wi + ui 2 — ui 3 + W4) 
(—wi + W2 — U13 + ui4)(— wi — W2 + wn + W4) 
(wi — W2 — W[i + W4)( — W\ — W2 — W3 + W4) . 



Thus, the determinant of the partition matrix does indeed factor into a brute-force 
iteration of every possible partition of W: the partition polynomial. 



2. Partitions and a System of Polynomial Equations 

The Partition problem determines if a given set of integers W — {w\, . . . ,w n } can 
be divided into two sets, / and W \ I such that ^2 we j w = ^2 we w\i w - I n this section, 
we describe a system of polynomial equations that models this question, and discuss the 
degree and monomials in an associated minimum-degree Nullstellensatz certificate. 

Proposition 1. Given a set of integers W = {wi, . . . ,w n }, the following system of 
polynomial equations 

xf — 1 = , for 1 <i <n , 

n 

w i x l = . 

has a solution if and only if there exists a partition ofW into two sets, I C W and W\I, 
such that J2 weI w = Y, we w\i w ■ 

Proof: The variables Xi can take on the values of ±1. Thus, we relate partitions to solu- 
tions by placing integers Wi with +1 Xi values on one side of the partition and integers 
Wi with —lxi values on the other. □ 
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Let [n] denote the set of integers {1, ...,n} and let denote the set of fc-subsets 
of [n]. For S € S%, let x s denote the corresponding square-free monomial of degree 
k in n variables. For example, given S = {1,3,4} C [5], the corresponding monomial 
— X1X3X4. Additionally, let S^ 1 denote the fc-subsets of [n] \ i . 



x 



Theorem 2. Given a set of non-partitionable integers W = {wi, . . . , w n } encoded as 
a system of polynomial equations according to Prop. 1, there exists a minimum- degree 
Nullstellensatz certificate for the non-existence of a partition of W as follows: 



1 = £( E E -i) + (£ E ^ s )( 



22 WiX, 



i—1 k even qn\i fc oddS^S^ i—1 

fc<n-l oto <= k<n 

Moreover, every Nullstellensatz certificate for the system of equations defined by Prop. 1 
contains one monomial for each of the odd parity subsets of each S%, and one monomial 
for each of the even parity subsets of each Su . 

Via Thm. 2, we see that the degree of the certificate is n for n odd, and n — 1 for n 
even. Furthermore, by considering the monomials present in the certificate as identifying 
the integers present on one side of a partition, we see that the monomials represent a 
brute- force iteration over every possible partition of W . We note that we identify the 
constant terms Cj with the case of placing every integer on one side of the partition and 
the empty set on the other. Thus, this result is similar to the independent set result (De 
Loera et al., ([8])) reviewed in the introduction. However, in this paper, we are interested 
not only in a combinatorial interpretation of the monomials, but also in a combinatorial 
interpretation of the unknowns c^s, 65. 

The proof of Thm. 2 is virtually identical to the proof of the independent set result 
described in (jSJ. Thus, we omit the formal proof here and state only the result: the 
Nullstellensatz certificates associated with this simple formulation of the NP-complete 
Partition problem are both hard to find and hard to verify. 

Example 1. The set of integers W = {1,3,5,2} is not partitionable. We encode this 
problem as a system of polynomial equations as follows: 

x\ - 1 = , x\ - 1 = , x\ - 1 = , x\ - 1 = , x x + 3x 2 + 5.T3 + 2x A = . 

Since W is not partitionable, this system of equations has no solution, and a Nullstel- 
lensatz certificate exists. Here is the minimum- degree certificate described by Thm. 2: 



155 842 188 908 \ , 2 . / 1 842 188 

1 X2X3 X2X4 -I £'3X4 (Xi — 1) + 1 X1X3 X1X4 

693 3465 693 3465 J y 1 ; \ 231 1155 231 



•ix^j{xl -1)+ ( - 

292 \, 2 s / 467 842 908 292 \, 2 , / 68 376 

H X3X4 (x.) — 1) + 1 xi £2 H X1X4 H £9x4 fx--! — 1) + xi X'2 

1155 V 2 ' V 693 693 693 693 V ' V 693 693 

1816 584 \, 2 , /155 1 467 34 842 

H xi X3 H X2X3 tx A — 1) + xi H X2 H X3 H X4 X1X2X3 

3465 3465 J K 4 ; V693 693 3465 693 3465 

188 908 292 

+ 693 

Note that the coefficient for (x\ — 1) contains only even degree monomials that do not 
contain x\ (similarly for [x\ — 1), etc. ) and that the coefficient for [x\ + 3^2 + 5x3 + 
2x4) contains every possible odd degree monomial in four variables. The combinatorial 
interpretation of a number such as 34/693 is explicitly demonstrated in Ex. 9. □ 



908 292 \ . 

xix 2 x 4 - — — X1X3X4 - — — x 2 x 3 x 4 (xi + 3x 2 + 5x 3 + 2x 4 ) 
3465 3465 / 
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3. The Partition Matrix: Definition and Properties 



In this section, we explore the well-known connection between Hilbert's Nullstellensatz 
and linear algebra, in terms of the minimum-degree certificate defined in Thm. 2: 

n n 

! = E( £ £ 1) + (£ £«* 



k<n-l ' 



k odd StzS^ 



We begin by defining graded reverse lexicographic order. We then construct a 2™ _1 x 2™~ 1 
square system of linear equations containing only the unknowns b. When ordered accord- 
ing to graded reverse lexicographic order, this square matrix is known as the partition 
matrix. We conclude by demonstrating that the partition matrix is not only symmetric, 
but also has a variety of properties essential to proving our main result in Section 4. 

3.1. Graded Reverse Lexicographic Order as a Tree 

Since we are dealing only with square-free monomials, we define graded reverse lexi- 
cographic order (denoted }zd) as follows. Given S € S%, we represent S as a vector in 
{0, 1}™ (denoted vs) by setting vs[i] = 1 if i € S and vs[i] = otherwise. For example, 
let S = {2, 3, 7} e Si Then v s = {0, 1, 1, 0, 0, 0, 1}. Given distinct S e S% and 5" e S£,, 
then S (Zd S' in two cases: 1) if k > k' , or 2) if fc = fc' and the right-most nonzero entry 
of vs — v s < is negative. For example, {2, 3, 4, 5} >d {1, 2, 5}, and {2, 3} hr> {1, 4}. 

In order to prove specific properties of the partition matrix, we use a slightly less com- 
mon, recursive definition of graded reverse lexicographic order. First, we order the ( "jj 
subsets of [n] in lexicographic order, creating sets Si, . . . , S n . Next, the sets Si, . . . ,S n 
are iterated, and for each Si, the ("I^) subsets of Si are iterated in lexicographic order, 
etc.. This order is pictorially represented as a tree in Ex. 2. 

Example 2. Here we pictorially order the set of integers [5] according to ^d- 

12345 



1234 




1235 1245 1345 

135 235 145 245 345 

I I 

35 45 



2345 



Using this tree, if two sets S, S' are from different levels in the tree with S higher 
than S' , then S hr> S' . For example, {1245} {234}. Additionally, if S,S' are from 
the same level in the tree but S appears further to the left than S' , then S }zd S' . For 
example, {23} >d {15}- Additionally, observe that if the even and odd cardinality subsets 
of [5] are iterated in> D order, then the following pairing of even and odd subsets occurs: 



12345 


123 


124 


134 


234 


125 


135 


235 


145 


245 


345 


1 


2 


3 


4 


5 


1234 


1235 


1245 


1345 


2345 


12 


13 


23 


14 


24 


34 


15 


25 


35 


45 






Given a set S in the pairing diagram above, i/5 £ 5, then S is paired with S\5. Ifb^S, 
then S is paired with S U 5. This observation is proven in general in Prop. 3.3. □ 



■5 



We refer to this tree as the order tree of [n] . If two sets 5*, 5" are children of the same 
parent in the tree, we say that the sets are contained in the same block. For example, 
{1, 2, 3} and {2, 3, 4} are in the same block, but {2, 3, 4} and {1, 2, 5} arc not. 

3.2. The Partition Matrix 

In this section, we demonstrate how to extract a 2" _1 x2 n_1 matrix from the minimum- 
degree certificate of Thm. 2. We begin by considering the coefficients of [x\ — 1): 

( E E 0- 

k even ci£- c n \i 

k<n—l * 

We observe that each monomial Ci^sx s multiplies (x\ — 1), which implies that each c^s 
appears in two equations (one corresponding to the monomial x s x\ , and one correspond- 
ing to the monomial — x s ). Thus, the unknown c^s appears in the first equation with a 
positive coefficient, and the second equation with a negative coefficient. This allows us to 
sum the two equations, and cancel the c unknowns in a cascading manner. For example, 
there is always one equation for the constant term: 

_c l,0 _ c 2,0 - ■ ■ ■ - c n,0 = 1 • 

Notice that this equation sums to one, since the Nullstellensatz certificate simplifies to 
one. There is also always one equation for each x\ monomial: 

hwi + c it $ = . (1) 

The biWi term appears in these equations since the product of 

(EE ^ slS ) ( y w i x i H \-w n x n ^ , 

k odd S£S^ 

k<n 

contributes the term 6^ • = biWiX^, among others. Notice that Eq. 1 sums to 

zero, since every monomial other than the constant term must cancel in a Nullstellensatz 
certificate. This set of n + 1 equations yields the following subsystem: 

— cu — c 2 j — • • • — c„ j = 1 , (constant term) 
b 1 w 1 + ci,0 = , (x\) 

b n w n + c„ i0 = . {x 2 n ) 
Summing these n + 1 equations together yields the following equation (in b only) : 

n 

E b ^ W i = 1 ' 
i=l 

In general, let SC. [n]\ibe an even cardinality subset, and consider the two monomials 
x s x\ and x s . Then, the following n— l^l + 1 equations are always present in the extracted 
linear system: 

bsui w i + c i.s — , { xS ' x i ) j f° r eacn i ^ S (2) 
E b s\i w i ~ E c *> s = ' • 
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Summing up these n— \S\ + 1 equations together yields the following equation (in b only): 

bsujWj + b s\jWj = • 
its jes 

Definition 1. Given a set of integers W — {wi, . . . ,w n }, the coefficient matrix of the 
following square system of linear equations 

^ bsujWj + ^ b s \jWj — , for each S G (S^ \ 0) with \S\ even 
its jes 

n 

Y h * W i = 1 ' 
i=l 

defines a 2™~ 1 x 2" _1 matrix with columns indexed by the unknowns bs (corresponding to 
the 2 n ~ 1 odd cardinality subsets of [n]), and rows indexed by the sets S (corresponding to 
the 2 n ~ 1 even cardinality subsets of[n], including $). This matrix is the partition matrix, 
denoted by Part n , with rows and columns ordered by graded reverse lexicographic order. 

By studying Eq. 2, we see that each c unknown appears in exactly one equation along 
with exactly one b unknown. Thus, solving for the b unknowns uniquely determines the 
entire certificate, and determining whether or not a given set W is partitionable depends 
entirely on the determinant of the partition matrix. 

Example 3. Let W = {wi,w 2 ,w 3 }. Via Thm. 2, the Nullstellensatz certificate is: 

1 = (ci,0 + Ci^ 23 }X 2 x 3 ){x\ - 1) + (c 2 ,0 + c 2 .{ 13 }X 1 x 3 )(xl - 1) + (c 3i + c 3 ^ {12 }X l x 2 ){xl - 1) 
+ (hxi + b 2 x 2 + b 3 x 3 + bi 23 xix 2 x 3 )(wixi + w 2 x 2 + w 3 x 3 ) . 

If W is not partitionable, there must exist an assignment to the unknowns c and b such 
that the certificate simplifies to one. In other words, the following system of linear equa- 
tions has a solution: 



(*?) 




c l,0 - 




= 


(x 2 x 3 ) 


~ c l,{23} + b 2W 3 + b 3 w 2 


= 






c 2,0 " 


- b2W2 


= 


(xiX2xl) 


c 3,{12} + bl2 3 W 3 


= 


(*!) 




c 3,0 " 


\- b 3 w 3 


= 


{xix\x 3 ) 


c 2,{13} + bl 23 W 2 


= 


(X1X2) 


~ c 3,{12} " 


- blW2 - 


- b2U>l 


= 


(xjx 2 X 3 ) 


c l,{23} + bl2 3 Wi 


= 


{x\X 3 ) 


~ C 2,{13} - 


h biW3 ~ 


h b 3 w 1 


= 


{constant term) 


— Cl,0 - C 2 ,0 - C 3 


= 1 



Following the simplifications described above, we extract a square system of linear equa- 
tions that contain only the b unknowns from these equations: 

bi 23 w 3 + bxw 2 + b 2 wx = , S = {1,2}, h 23 w 2 + bxw 3 + b 3 w x = , 5" = {1,3}, 
b\ 23 wi + b 2 w 3 + b 3 w 2 = , S = {2, 3} , h w\ + b 2 w 2 + b 3 w 3 = 1 . 5" = . 

Ordering the columns as {&123, b\, b 2 , b 3 }, the partition matrix is as follows: 

6123 h b 2 b 3 

{1, 2} w 3 w 2 w\ 
{1, 3} w 2 w 3 Wi 
{2,3} wi w 3 w 2 

I»i I»2 w 3 
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As a preview of our main result, we note that the determinant of this matrix is 

(ioi + w 2 + w 3 )(-wi + w 2 + w 3 )(wi - w 2 + w 3 )(-wi -w 2 + w 3 ) , 

which represents a brute-force iteration over all of the possible partitions of the set W . 
This will be formally defined as the partition polynomial in Sec. 4- n 

For the duration of this section, we collect and prove a variety of interesting facts 
about the partition matrix Part„. The properties may seem particularly intricate and 
overly specific, but each property is necessary in the proof the main result in Sec. 4. 
Throughout these propositions, the notation Part„[i, :] follows standard Matlab nota- 
tion, and denotes the i-th row, whereas Part„[i, 1 : j] denotes the i-th row with columns 
1 through j only. 

Example 4. Here we display the 16 x 16 partition matrix Parts . 





12345 


123 


124 134 234 


125 


135 235 


145 245 345 


1 


2 


3 


4 


5 


1234 


w 5 


W4 


w 3 


102 


Wl 



































1235 


W 4 


W 5 











w 3 


10 2 


Wl 


























1245 


W 3 





w 5 








W4 








W2 


Wl 




















1345 


»2 








10 5 








W4 





w 3 





Wl 

















2345 


Wl 











11I5 








W4 





w 3 


W2 

















12 





w 3 


tt>4 








w 5 

















w 2 


Wl 











13 





W2 





104 








w 5 














w 3 





Wl 








23 





Wl 








W4 








W5 














w 3 


W2 








14 








W2 


10 3 














W5 








W4 








Wl 





24 








Wl 





W3 














W5 








W4 





W2 





34 











Wl 


102 

















w 5 








W4 


w 3 





15 

















102 


103 





W4 








w 5 











Wl 


25 

















Wl 





w 3 





W4 








w 5 








W2 


35 




















Wl 


W2 








W4 








w 5 





w 3 


45 


























Wl 


W2 


W3 











W5 


W4 






































Wl 


W2 


W3 


W4 


w 5 



Proposition 3. The 2™ _1 x 2 n_1 matrix Part n has the following properties: 

(1) The entry Wi with i = {1, . . . , n} appears exactly once in each row and column. 

(2) The submatrix Part n [l : n, 1 : n] consists of w n on the diagonal, the entries 
w ni . . . ,w\ in both Part n [l : n, 1] and Part n [l, 1 : n], and zero elsewhere. 

(3) If row i is indexed by set S C [n] (with \S\ even), and n <G S , then column i is 
indexed by S\n. If n ^ S, then column i is indexed by S Lin. 

(4) All diagonal entries of Part n are equal to w n . 

(5) Given any row (column), the entries to the left (above) the diagonal are indexed by 
SLij, and the entries to the right (below) the diagonal are indexed by S\j. 

(6) Part n is symmetric. 

Proof of Prop. 3.1: After inspecting the equation defining the partition matrix 

bsujwj + Y bsxjWj = , 
its jes 

where S represents an even cardinality subset of [n], it is evident that each row contains 
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exactly one entry for each Wi. To see that each column also contains exactly one entry 
Wi with i = {l,...,n}, consider the column indexed by unknown bs> where S' C [n] 
with odd cardinality. Then, for each j £ S' , the row indexed by S = S' \ j contains Wj. 
Additionally, for each j £ S' , the row indexed by S = S' Uj also contains Wj. Thus, each 
row and column contains exactly one entry Wi for i = {1, . . . , n}. □ 

Proof of Prop. 3.2: We will prove that Part„[l : n, 1 : n] has the following form: 



w„ 


Wn-1 W n -2 




Wl 


Wn-1 


W n 







Wn-2 


w„ 








. Wl 








w n 



For n even, row 1 is indexed by {1 • • • n} and the first n columns are indexed by the ( "j) 
subsets of [n] in lexicographic order. Moreover, rows 2 though n are indexed by the ("I^) 
subsets of {1 • • • (n — 1)}. Therefore, the claim holds by inspecting the equation defining 
the partition matrix (Def. 1). For n odd, column 1 is indexed by {1 • • • n} and the first n 
rows are similarly labled with the ( "J subsets of [n]. Moreover, columns 2 through n 
are indexed by the ("Ig) subsets of {1 • • • (n — 1)}. □ 
As an example of Prop. 3.3, note that row {12} is paired with column {125}, and 
column {1} is paired with row {15} in Ex. 4. 

Proof of Prop. 3.3: To prove this claim, suppose that we have the "order tree" T„_i for 
the subsets of [n — 1] . In order to create the order tree T n for the subsets of [n] , we first 
copy T„__i and add the integer n to each set, creating the tree T n _i U n. We then join 
the node in T„_i U n indexed by {1 ■ ■ ■ n} to the node in T„_i indexed {1 • • • (n — 1)}. 
The resulting tree is the order tree for T n . For example, 



1234 



12345 




1255 1 245 1345 2 3 45 

/\ I 
15 235 145 245 345 
I I 
3i 45 



Since no set in T„_i contains the integer n and every set in T„_i U n contains n, it 
is easy to see that the even and odd sets are paired by inspecting how T„_i overlays on 
top of T„_i U n. Thus, the claim holds. □ 
Proof of Prop. 3.4 ■ This result follows from the equations defining the partition matrix, 
and also Prop. 3.3, which defines the row-column pairing of the diagonal element. □ 
Proof of Prop. 3.5: Consider a row indexed by the set S. Since SUj >d S\j, we can be 
certain that the entries formed by S U j appear to the left (above) the entries formed by 
S\j. The diagonal is either formed by S U n (if n ^ S), or S \ n (if n € S). We observe 
that S U j h D S U n (if n £ S), and S \ n S (if n e S). □ 
Proof of Prop. 3.6: Consider an arbitrary row i indexed by a set Si, and let column i be 
indexed by the set bj. In order to prove symmetry, we must show that row i is equal to 
column i. By Prop. 3.3, n is either in Si or hi, but not both. Without loss of generality, 
assume n € Si and bi = Si\n (e.g. Si = {15} and bi = {1}). Suppose Part„[z, j] = Wfc j for 
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j < i. We must show that Part„[j, i] = w kj . Since j < i, kj ^ Si, and column j is indexed 
by bj = Si U kj (e.g. in row {15}, W2 appears in column {125}). Since n € {Si U kj), row j 
is indexed by Sj = (SiUkj)\n (e.g. Sj = {12}). Then, Sj\kj = bi, and Part„[j, i] — Wkj ■ 

Suppose i < j, and Part n [i,j] is again equal to w^- Then kj € Si, and column j is 
indexed by bj — Si \ kj and row j is indexed by Sj — (Si \ kj) \ n (e.g., in row {15}, wi 
appears in column {1}). But then Sj U kj — bi, and Part„[j, i\ = Wk r 

A similar argument holds if n (ji Si, but with the logic reversed. Since we have shown 
that Part„[i, j] = Part„[j, i\. we have shown that the matrix is symmetric. □ 

Proposition 4. Let column j be indexed by set {ji ■ ■ ■ jk} 7 where j\, . . . , j r are consec- 
utive ascending integers excluding n, and let column j + 1 < 2™ _1 be indexed by set 
. . . , j' k ,} where j[, . . . , j' r , are consecutive ascending integers. 

(1) If ji = 1, then 3 a row i such that Part n [i : (i + r — = {w r , . . . , u>i}. 

(2) If ji = j'i = 1, then 3 a row i such that Part n [(i + r) : (i + 2r — 2), j + 1] = 

{w r -!,...,Wi} . 

(3) If ji 7^ 1, then let 1, ... ,t be the integers such that t = ji — 1. Then, 3 a row i' 
such that Part n [i' , (j — t) : j] = {wi, . . . , Wt+i}. 

(4) If ji ^ 1 and j[ = 1, then, given the rows i and i' that exist by claims (1) and (3), 
respectively, i = i' + 1. 

Before proving claim 4.1, we study an example. In Ex. 4, column 1 is indexed by 
set {12345}, and thus contains a set of consecutive ascending integers (excluding n) 
equal to {1,2,3,4}. Notice that Parts[2 : 5, 1] = {104, W3, W2, wi}- Additionally, column 
2 is indexed by set {123}, and thus has a consecutive ascending set of integers equal 
to {1,2,3}. Furthermore, Part 5 [6 : 8,2] = {w 3 , w 2 , u>i}. Thus, for column 1, the row i 
referred to in claim (1) is i = 2, and for column 2, the row i referred to in claim (1) is i = 6. 
Given a column j, we will refer to the entries in Part n [i : (i + r — = {w r , . . . , wi} 
as the "trailing block" of column j. 

Example 5. Here we highlight the trailing blocks of each applicable column. 
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Proof of Prop, ^.l: If the set {ji • • • jk} indexing column j has j\ — 1, then the set of 
consecutive ascending integers {ji ■ ■ ■ j r } is {1 • • • r}. Recall that odd order sets index the 
columns, and even order sets index the rows. Thus, the even order sets 

{!■ ■ - rjr+i - ■ ■ jk}\r , 

{1 ■ ••rjV+i • ■ -j k } \ (r - 1) , 

{1 • --rjr+i ■■■jk}\l , 

index the consecutive block of rows from some row i to row i + r — 1. By the definition 
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of the partition matrix Part„ (Def. 1), if column j is indexed by set 

{l---rj r+1 ---j k } 

and row i is indexed by set 

{l---rj r+1 ---j k }\r , 

then Part„[i, j] — w r , since r is contained in the column set, but not contained in the 
row set. Thus, Part„[i : (i + r — = {w r , . . . , wi} and the claim is proven. □ 

As an example of claim 4.2, we recall in Ex. 4, that column 1 is indexed by set {12345}, 
containing consecutive ascending integers (excluding n) equal to {1,2,3,4}. Therefore, 
r = 4, since there are four consecutive ascending integers. Recall further that column 
2 is indexed by set {123}, containing a consecutive ascending set of integers {1,2,3}. 
Therefore, r' = 3, since there are three consecutive ascending integers within the label. 
Observe that the size of the consecutive ascending set of integers in column 2 is one less 
then the size of the consecutive ascending set of integers in column 1 (excluding n). We 
claim in Prop. 4.2 that Part„[(« + r) : (i + 2r — 2), j + 1] = {uv-i, and we observe 

that Part 5 [(2 + 4) : (2 + 2 • 4 - 2), 1 + 1] = Part 5 [6 : 8,2] = {w 3 ,w 2 ,wi}. Thus, claim 
(4.2) holds in this example. 

Proof of Prop. By Prop. 3.1, Part„[i : (i + r — = {w r , . . . ,w\} . Assume 

that k = k' . In this case, since the set indexing column j and the set indexing column 
(j + 1) both contain 1, both sets are contained within the same block (have the same 
parent within the order tree). Therefore, by the lexicographic ordering present within a 
given block, column j is indexed by {1 • • • rj r+ i ■ ■ • j^} and column j + 1 is indexed by 
{j[ ■ ■■f r >3 / r , +1 ■■■&} = {!■■■ (r-l)(r + l)f r+1 ■ ■ ■ f k ,} = {1 • • • (r - l)(r + l)j r+ i • ■ • j k }. 
For example, suppose column j is indexed {1236789}. Then column j + 1 must be indexed 
by {1246789} according to lexicographic order. In other words, r' = r — l,j' r , +1 = r + 1 
and j' r , +2 = jr+i- Therefore, the even order sets 

{l...(r-l)(r + l)j r+1 ---j fe }\(r-l) , 
{l-"(r-l)(r + l)j r+ i-"i fc }\(»--2) , 

{l---(r-l)(r + l)j r+ i-"ifc}\l . 

index the consecutive block of rows from row i+r to row i+r+(r— 1) — 1. Therefore, by the 
definition of the partition matrix (Def. 1), Part„[(i + r) : (i + 2r — 2), j] = {w r -i, . . . ,w±} 
and the claim is proven. 

If fc' + fc and ji = j[ = 1, then this is the special case that only occurs in column 1 
and column 2 when n is odd. In this case, column 1 is indexed by {1 • • • n), and column 
2 is indexed by {1 ■ ■ ■ (n — 2)}. Additionally, the first n rows are indexed by the („™ 1 ) 
subsets, and the next n — 1 rows arc indexed by the (™I 2 ) subsets of [n— 1]. Therefore, 
Part„[2 : n, 1] = {w n -\, . . . , w\] and Part„[(n + 1) : (2n - 2), 2] = {tu n _2, • • • , w\} as 
claimed. In either case, the claim holds. □ 

As as example of Prop. 4.3, we observe that in Ex. 4, row {45} contains entries 
{wi,w 2 , w 3 } in the columns indexed {145}, {245}, {345}, respectively. 

Proof of Prop. 4-3: Since column j is indexed by {ji • • • jk} and ji ^ 1, this means that 
the last row of column j containing an entry is the row indexed by the set {ji • • • jk}\ ji = 
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{ji • • • jk}- Furthermore, the columns indexed by 

{O'l - l)j2 ■ ■ ■ Jfc} , 

{(ji-2)j2---ifc} , 

{l.?2"-Jfc} , 

are the columns j — l,j — 2, . . . , j— ji + 1, respectively. Thus, the row indexed by {j 2 •••jk} 
contains entries w\,...,Wj 1 in columns j — j\ + 1, ... ,j — l,j, respectively, and the claim 
is proved. □ 

In order to clarify the statement of the Prop. 4.4, we observe that in Ex. 4, the col- 
umn labled with set {345} has the row indexed by set {45} (i' = 15), and the next 
column (indexed by {1}) has the row indexed {0} (i = 16). Therefore, i = i' + 1 and 
|{345}| - |{1}| = 3 - 1 = 2, and the condition holds. 

Proof of Prop. In this claim, j\ ^ landjj = 1. We begin by assuming that k' = fc— 2. 
In this case, column j and column j + 1 are indexed by sets from two different levels of the 
order tree. Furthermore, since j is indexed by the set at the end of a level, n = jk, and 
{ji ' ' 'jk} is actually equal to {(n — k + 1) ■ ■ ■ n}. This implies that the row i' from claim 3 
is indexed by {ji ■ ■ ■ jk}\ji, or {(n — k+1) ■ ■ ■ n}\(n— k+1). Since the set labeling column 
j + 1 is the first set in the next odd level, {j[ ■ ■ ■ j' k ,} = {1 • • • fc'}. Therefore, the row i 
from claim 1 is {1 • • • k'} \ k' = {1 ■ ■ ■ {k - 2)} \ (k - 2). But clearly {1 • • • (k - 2)} \ (k - 2) 
labels the row directly after {(n — k + 1) ■ ■ ■ n} \ (n — k + 1). Therefore, i = i' + 1 as 
claimed. 

In the second case, we assume fc' = fc. Thus, the two sets labeling columns j and j + 1 
respectively are in the same level of the order tree. However, since j\ ^ 1 and j[ = 1, the 
two sets are from different blocks in the same level. Moreover, the set labeling column j 
is at the end of one block, and the set labeling column j + 1 is at the beginning of the 
next block. Therefore, the set labeling j is {ji ■ ■ ■ j r jr+ijk}, and the set labeling column 
j + 1 is {1 • • • (r — l)(j r + l)jr+i • • • jk}- Thus, we see that the row i! from 4.3 associated 
with column j is indexed by {ji • • • jrjr+ijk} \ ji- The row i from 4.1 associated with 
column j + 1 is indexed by {1 • • • (r — l)(j r + l)j r +i • • • jk} \ ( r — !)• These two rows, 

{ji ■ ■■ jrjr+ijk} \ji , and 
{l---(r-l)OV + l)jV+i"-jfc}\(r-l) 

are lexicographically adjacent. Therefore, i = i' + 1 as claimed. 

Since the claim holds in both cases, 4.4 is proved. □ 

The following three propositions deal with identifying similar blocks of entries in dif- 
ferent rows and columns. Within these propositions, we will use the following notation. 
Given a column j, let Aj denote the set of entries above the diagonal. Notice that the 
set Aj C {tui, . . . , w n _i} since w n is never in Aj because w n is always on the diagonal. 

Proposition 5. Let j be a column indexed by {lj 2 ■ ■ ■ jk}- Let i be any row in the trailing 
block of column j. Then, the non-zero entries in Part n [i, {j + 1) : (i — 1)] are the same 
as the entries in Aj. 
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Example 6. Here is an example of Prop 5. Let j be column {124}. Then the trailing 
block of column j consists of entries {w2, u>i}, and Aj — {w 3 }. Let i be either row in the 
trailing block, either {14} or {24}. Then the only non-zero entry in Part n [i, (j + 1) : (i— 1)] 
is {w 3 }, which is exactly equal to Aj. 
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□ 



Proof of Prop. 5 : Given a column j indexed by {I72 ■ ■ ■ jk}, let i be a row in the trailing 
block of column j such that P[i, j] — w t . Then row i is indexed by set {lj 2 • • • jk} \ t- We 
will show that the set of non-zero entries in Part n [i, (j + 1) : (i — 1)] are the same as the 
set of entries above the diagonal in column j, denoted by Aj. 

We will prove both directions of the inclusion. First, we will prove that Aj C Part„[i, (j+ 
1) : (i — 1)]. Let w a € Aj. Then a ^ { 1 J2 • ■ • Jfe } ■ Moreover, the column indexed by 
{{1.72 • • ■ jk} \t}Un falls between columns j and the diagonal. Therefore, the entry 
Part n [i, i'] = w a where j + 1 < i' < i — 1. 

Conversely, let Part„ [i, i'] = w a where j + 1 < i' < i - 1. Then, a £ {lj 2 ■ ■ ■ jk}, which 
implies w a e Aj. 

Since we have proven both directions, the non-zero entries in Part„[i, (j + 1) : (i — 1)] 
are the same as the set Aj. □ 

Proposition 6. Let j be a column indexed by {I72 ■ ■ ■ jk}- Let i be any row in the trailing 
block of column j. Let j < i' < i be such that P[i,i'] is non-zero. Then Aj U P[i,j] = 

^uP[i,i'] . 

Example 7. Here is an example of Prop 6. Let j be column {124}. Then the trailing 
block of column j consists of entries {w 2 ,Wi}, and Aj = {w 3 }. Let i be row {14} in the 
trailing block, and let i' be column {234} . Then AjUP[i,j] — {w 3 }L){wi} — {wi}U{w 3 } = 
A v UP[i,i'\. 
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□ 



Proof of Prop. 6: Let column j be indexed by {lj 2 ■ ■ ■ jk}, and let i be a row in the trailing 
block of column j such that P[i, j] = Wt- Then, row i is indexed by set {lj 2 • • • jk} \ t. If 
row i contains a non-zero entry P[i, i'] — w a before the diagonal, then column i' is indexed 
by {{lj'2 • • • jk} \t}llo. Since a is not in {lj 2 ■ ■ ■ jk}, then a is in Aj. Moreover, since t is 
not in the set labeling the column i' , then t £ A^. Since every other index between the 
two column labels ({I72 ■ ■ ■ jk} and {{I72 ■ ■ ■ jk}\ t} U a) is the same, Aj Uw t = Ay U w a , 
or A,- U P[i, j] = A v U P[i, □ 

Proposition 7. Le£ j be a column indexed by {lj 2 ■■■jk}- Let i be any row in the trailing 
block of column j. Let i < i' . Then Aj U P[i, j] U P[i, i'] = Ai> . 

Example 8. Here is an example of Prop 7. Let j be column {124}. Then the trailing 
block of column j consists of entries {u>2,u>i} 7 and Aj = {W3}. Let i be row {24} in the 
trailing block, and let i' be column{A}. Then, AjUP[i,j]UP[i,i'] = {u>3}U{u>i}U{u>2} = 
{101,102,103} = A v . 
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Proof of Prop. 7: Let column j be indexed by {I72 ■ ■ ■ jk}, and let i be a row in the trailing 
block of column j such that P[i, j] = Wt - Then, row i is indexed by set {lj 2 ■ • ■ jk}\t- If row 
i contains a non-empty entry P[i, i'] — w a after the diagonal, then column %' is indexed 
by {{I72 • • • jk}\t} \o. Since a is in {lj 2 ■ ■ ■ jk}, then a is not in Aj, but it is in A41 (since 
it is above the diagonal in column i'). Moreover, since t is not in the set labeling column 
i', then t e Ai>. Since every other index between the two column labels ({I72 ■ ■ ■ jk} and 
{{lj 2 • • - jk} \ t}\a) is the same, Aj U w t U w a — A v , or Aj U P[i,j] U P[i, i'\ — A v . □ 
Having gathered together a series of facts about the partition matrix, we will now 
investigate the determinant of the partition matrix. 



4. The Partition Matrix and Partition Polynomial 

Given a square non-singular matrix A, Cramer's rule states that Ax = b can be solved 

by 

det(A\j) 
1 det(A) ' 

where ^4 1 ^ is the matrix A with the z-th column replaced with the right-hand side vector 
6. In Section 3, we extracted a 2 n_1 x 2™ _1 square linear system from the general linear 
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system constructed via the minimum-degree Nullstellensatz certificate described by Thm. 
2. Here, we see by Cramer's rule that the unknowns within that certificate are ratios 
of two determinants. In this section, we show that the determinant of the partition 
matrix is equivalent to a brute-force iteration over all the partitions of W. Therefore, the 
denominator of any unknown in the certificate is a combinatorial representation of the 
partition problem. 

We observe that, in general, the linear system Ax = b may have a solution even if 
det(A) — 0. However, in the case of the partition matrix, when we demonstrate that the 
det(A) is equal to the partition polynomial, we will be demonstrating that Ax — b only 
has a solution in the case when det(A) ^ 0. 

Let {0, 1}™ be the set of all 0/1 bit strings of length n. For S G {0, 1}™, let Si denote 
the z-th bit in the string S. 

Definition 2. Given a set W — {wi, . . . , w n }, let 



se 

be the partition polynomial ofW. 

For example, let n = 5, and S G {0, l} 4 be S — "1011". Then, S corresponds to the 
— wi + W2 — W3 — Wi, and denotes a partition of W = {wi, . . . , w§}, with w§ fixed on the 





+ 






W3 




W4 





If this arrangement of Wi is a partition of W, then — wi + w 2 — w 3 — W4+W5 = 0. In 
this way, any bitstring S G {0, l} n_1 is equivalent to fixing w n on the "positive" side of 
the partition, and then arranging the other Wi on the "positive/negative" side, according 
to sign. In this way, the partition polynomial represents an iteration over every possible 
partition of W, avoiding double-counting by permanently fixing w n on the "positive" 
side. If the set W is partitionable, one of bitstrings S will define a factor of the partition 
polynomial that sums to zero. We will show that the determinant of the partition matrix 
is the partition polynomial: therefore, if the determinant of the partition matrix is zero, 
the linear system has no solution, and there is no Nullstellensatz certificate. 

Example 9. In Ex. 1, we presented an actual minimum- degree certificate for the non- 
partitionable set W = {1,3,5,2} . We observe that 

-51975 = (1 + 3 + 5 + 2)(-l + 3 + 5 + 2)(1 - 3 + 5 + 2)(1 + 3-5 + 2) 

(-1 - 3 + 5 + 2)(-l + 3 - 5 + 2)(1 - 3 - 5 + 2)(-l - 3 - 5 + 2) . 

Via Cramer's rule, we see that the unknown 64 is equal to 

-2550 _ 34 
4 ~ -51975 ~ 693 ' 

which is indeed the value of unknown 64 as it appears in the certificate. □ 
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Example 10. Here is the determinant of the 8x8 partition matrix Part^: 



(\ 




det 



VL 



(1)3 1«4 W2 W\ 

W2 W4 W3 wi 

i»i W4 w s W2 

w 2 w 3 W4 t»i 

wi W3 W4 W2 

wi W2 1114103 

l«i ID 2 !«3 1»4 



(wi + U>2 + W S + W4,)(—W\ + U>2 + W3 + «J 4 ) 
(till - W2 + W S + M>4)(«>1 + Ul 2 — UI3 + U) 4 ) 
(— t«l + W2 — Wz + W4j{—W\ — W2 + W3 + W4) 
(wi — U>2 — 1«3 + U) 4 )( — Uli — W2 — W3 + W4) . 



Theorem 8. Given W = {w\, . . . , w n }, the determinant of the partition matrix of W is 
the partition polynomial ofW. 

We will prove this theorem in a slightly unconventional way. Recall that 



for annxn matrix A. Via this formula, we can see that the determinant of the partition 
matrix is a degree 2 n_1 polynomial. Therefore, if we can show that each of the 2™ _1 factors 
of the partition polynomial are also factors of the determinant of the partition matrix, 
we will have provided a unique factorization for the determinant, and thus illustrated 
that the determinant of the partition matrix is equal to the partition polynomial. 
Given S G {0, l} n_1 , if we make the substitution 



into the partition matrix, and then demonstrate that the determinant of the partition 
matrix after the substitution is zero, we will have shown that the factor defined by the 
bitstring S is a root of the determinant. If we can prove this claim for any bitstring 
S G {0, l} n_1 , we will have proven Thm. 8. Consider the following algorithm: 

ALGORITHM: FactorCheck 

INPUT: An integer n and a factor S G {0, l}™" 1 . 
OUTPUT: The reduced 2™- 1 x 2 n - 1 matrix P with factor S substituted for w n . 

I Set P = Part„ (Def. 1), set w n = - ( Yli=i(- 1 ) S>w i) in P > and set i = 2 ■ 

4 for j = 1 to 2™- 1 do 

5 while P(l,j) ^ do 

6 Set Wt — P , and set cur_sgn — — (sign of w t in cell P(l,j)) . 

8 Set P(l, :) = cur.sgn ■ P(i, :) + P(l, :) . 

9 Set i = i + 1 . 

10 end while 

II end for 

12 return P . 



n 



det(A) = sgn(a) J| A i>ff . . 



creS„ i=l 




1G 



For example, suppose P[l,j] = w\ + w 2 — w 3 , and P[i : (i + 2), j] = {w 3 , w 2 , w\}. 
Then, after applying the following row operations 

P[l,:]=P[i,:]+P[l,:] , 
P[l,:] = -P[i + 1,:]+P[l,:] , 
P[l,:] = -P[i + 2,:]+P[l,:] , 

the entry P[l, j] = 0. 

In Appendix A, we demonstrate an example of the algorithm running step-by-step 
on a particular integer n and input factor S. We will show that when this algorithm 
terminates, the first row of the output matrix P consists entirely of zeros. This would 
demonstrate that the factor described by the particular input bitstring S is a root of the 
determinant. Since we will show the result holds for any bitstring S, this will have shown 
that the determinant is the partition polynomial, and conclude the proof of Thm 8. 

When proving that an algorithm is "correct" (i.e., has the desired outcome), the 
technique is to define a loop invariant, and then show initialization, maintenance and 
termination: 

• Initialization: The invariant is true prior to the first iteration of the loop. 

• Maintenance: If the invariant is true before the current iteration of the loop, then 
the invariant is true before the next iteration of the loop. 

• Termination: When the loop terminates, the loop invariant provides a useful property 
that helps show that the algorithm is correct. 

Given a column j indexed by {j\ ■ ■ ■ jk}, let Tj = . . . , j r } C [n — 1] be the longest 
sequence of consecutive ascending integers (excluding n), or the trailing block, as in 
Prop. 4.1. Let Aj C [n — 1] be the set of entries above the diagonal in column j, and let 
Bj = [n — 1] \ (Aj U Tj) be the set of entries below the diagonal, excluding the trailing 
block. 

Example 11. In column {124}, we highlight the trailing block T124. In column {134} ; 
we highlight the set B134, and in column {1} we highlight the set A\. 
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Finally, let 

sgn(a, Aj) = I J (-l) 5c , for each a E Aj , 

e(EAj \a 

sgn(Aj) = [] (-l) 5c ■ 
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Note that sgn(a, A,-) and sgn(Aj) are both ±1. Since the FactorCheck algorithm 
contains the variables i,j,P and the input factor bitstring S, any loop invariant for the 
algorithm FactorCheck can be expressed in terms of those variables. We will also 
describe the conditions of the loop invariant in terms of the sets Aj , Bj and Tj . 

Lemma 1 (Line 4 Loop Invariant). On the j-th iteration of the for loop beginning on 
line 4, the matrix P satisfies the following properties: 

(1) If P[l, j] 7^ and column j is indexed by {ji ■ ■ ■ jk} where {ji ■ ■ ■ j r } is the longest 
consecutive set of ascending integers (excluding n). Then, 

P[l,j] = -(j2 s ^(A 3 )(-l) St w t ) , 

and P[i : (i + r — = {w r , . . . , w\} (the trailing block Tj for column j). 

(2) If P[l, j] — 0, then either 3 column j' such that P[l, 1 : (j' — 1)] = and condition 
(1) holds for column /, or i = (2 n_1 + l). 

(3) For column i , 

aeAi 

(4) For column i' such that i < i', let A' C A v . Then 

P[l,i'] = sgn(a,Ai>)w a . 

aeA' 

(5) For column i' such that j < i' < i, let B' C B v . Then 



P M = -( E sgn(A')(-l) S ^x ■ 

\a;e(T,/US') / 

(6) P[l, 1 : (j - 1)] = 0. 

These loop invariant conditions can be better understood by studying the example 
displayed in Appendix A, and seeing that each of these conditions hold during each iter- 
ation of the algorithm. 

Proof of Lemma 1 (Line 4 Loop Invariant): We must show that conditions (1) through 
(6) hold during initialization, and are maintained from one iteration to another. We must 
also show that when the FactorCheck algorithm terminates, it does so with "useful 
information", which in this case is conveyed by condition (6). 

Initialization: 



18 



By Prop. 2, the first n rows and columns of Part„ are always the following: 



n 



W n -1 W n _ 2 ■■■ W-i 

w n ••• 
w n : 



w 



w 



'n-2 



Wl 







•■ 

w, 



n 



Upon intialization of line 4, column j = 1, row i = 2, and the substitution 




has been performed in the matrix P. We will investigate each of the conditions of Lemma 
1, and show that they are true. 
(1) Upon initialization, P[l, 1] 7^ 0, and P[2 : n, 1] = {w n -\, . . . , w\}. Since there are 

no entries above the diagonal in column 1, sgn(Aj) = 1. Therefore, the trailing 

block Ti = {w n -i, . . . , wi}, and 



and condition (1) of Lemma 1 holds. 

(2) Since P[l, 1] 7^ 0, condition (2) does not apply. 

(3) For column i = 2, the entries above the diagonal arc A 2 = w n _\. Therefore, 
sgn(a,^4i) = 1 since A t \ = 0. Therefore, 



and condition (3) holds. 

(4) For columns i' = 3, . . . , n, the sets of entries above the diagonal are {wn-2}, ■ ■ ■ , {wi} 
respectively. Thus sgn(a, Ay) = 1 since A^i \w a = 0, as before. Furthermore, by 
Prop. 1, there is exactly one entry for each Wi in a given row. Therefore, the entries 
in row 1 for columns i' = n + 1, . . . , 2™ _1 are zero. Therefore, condition (4) holds. 

(5) There are no columns that satisfy j < i' < i when j = 1 and i = 2. 

(6) There are no columns less than j when j = 1. 

Thus, conditions 1 through 6 hold upon initialization. We will now show that each of the 
conditions are maintained from one iteration to the next. 

Maintenance: 

In order to prove maintenance, we assume conditions (1) through (6) hold at the j-th 
iteration, and we must show that conditions (1) through (6) hold at the j + l-th iteration. 
But whether or not the conditions hold depend soley on whether or not lines 6-9 in the 
inner while loop are executed. Thus, we break the proof of maintenance into two cases. 
The first case is when P[l, j) = (when lines 6-9 are not executed), and the second is 
when P[l,j] (when lines 6-9 are executed). 

Case 1: P[l, j] = 0. In this case, lines 6-9 are not executed. Thus, the matrix does not 



P[l,l] = -(£(-l) s * Wt ) , 



tar. 
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change between the j-th and (j + l)-th iteration. Furthermore, since P[l, 1 : (j — 1)] = 
(condition 6) and P = 0, columns P[l, 1 : (j + 1)] = 0. Finally, by Prop. 4.3, either 
i = 2 n+1 + 1, and we have already cancelled all the entries (see the second-to-last step 
of the example in Appendix A), or there is a column j' where the trailing block begins 
at row i. Thus, the conditions hold at the (j + l)-th iteration. 

Case 2: P[l,j] ^ 0. In this case, lines 6-9 are executed. Furthermore, by condition 
(1), P[i: (i + r- 1), j] = {w r , . . .,wi}. Since 

^[U] = -(£sgn(^)(-l)0 , 

the signs calculated in line 7 during the j-th iteration are 

(sgn^X-l) 5 "), (sgn(^)(-l) s ^) . . . , (sgn(^)(-l) s i) . 

Thus, after the r iterations of the while loop of line 5, the entries w r , . . . , W\ present in 
P[l,j] will be cancelled, and P[l,j] = 0. Therefore, P[i, 1 : j] = at the beginning of 
the (j + l)-th iteration, and condition (6) of the loop invariant holds. 

Now we consider an arbitrary row operation performed in line 8, and assume that 
conditions (3) through (5) hold before the operation, and prove that the conditions hold 
after the operation. Let i t be an arbitrary row with i < i t < (i + r— 1) and P[i tl j] — w t . 
Since we have already shown that the conditions hold on column j, we must only consider 
the effect of this operation on entries to the left of the diagonal (j < i' < i t ), the diagonal 
(i' = it), and the right of the diagonal (i t < i'). 

• (j < i' < i t ): Consider a non-zero entry P[i t ,i'] — Wb, with j < i' < i t (an entry to the 
left of the diagonal). By the loop invariant, 



P[l,i'} = -[ Yl sgn(^)(-l) £ "«'.. 



\xe(T 4 ,UB') / 

The row operation performed is sgn(Aj)(— l) St P[i t , i'} + P[l,z']. However, by Prop. 
6, Aj U P[it,j] = Aii U P[it,i']. Therefore, Wb G Aj and wt € Ay. In particu- 
lar, sgn{Aj){— l) St = sgn(Aii)(—l) Sb . Since the sign of Wb in P[l,«'] is equal to 
— (sgn(Ai')(— l) Sb ), this row operation cancels the entry Wb in P[l,i'], and condition 
(5) of the loop invariant holds. Moreover, we observe that, in the case where i' = j + 1, 
then every entry above the trailing block has been iterated, and therefore cancelled. 
Thus, either P[l, j + 1] = (if there is no trailing block), or every entry in Pj+i has 
been canceled and 

P[l,j + l] = -( E sgn(A, +1 )(-l)0 . 

Thus, we have shown that when j < i' < i t , all relevant conditions of the loop invariant 
hold. 

(i' = it): Now, we will consider the row operation performed on the diagonal, or 
sgn(Aj)(— l) St P[i t , it] + P[l,i*]. By condition (3) of the loop invariant, 

P[l,it] = s S^(a,A H )w a . 

aeA H 
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But, by Prop 5, we know that the entries in row it between columns j and it (excluding 
P[i t ,j]) are the same as the entries above the diagonal in column j. Furthermore, by 
the symmetry of the matrix (Prop 3.6), we know that P[it, 1 : it] — P[l ■ it, it]- Thus, 
the entries above the diagonal in column i t are the same as the entries P[i t ,j] U Aj. 
Thus, sgn(Aj)(— l) St = sgn(A it ). Recall also, that since P[i t ,it] is a diagonal element, 

P[iuH] = -(j2(-l) Si w^j 

Consider a particular w a € P[i t ,it] such that a g A it . Then, the sign of w a in P[i t ,i t ] 
is -((-1) 5 »). But 

-((-l) s °sgn(A it )) =-((-l)^(-l)^).sgn(a,A it ) = -sgn(a,A it ) . 

Thus, the row operation sgn(Aj)(—l) St P[i t ,it] + P[l,it] cancels all w a with a g Ai t , 
and 

P[l,it} = -( ]T sgn(A Jt )(-l) 5 ^J . 

\i£(T if UB, t ) / 

Finally, row i t is incremented to it + 1, and condition (5) of the loop invariant holds. 
Thus, when i' = it, all relevant conditions of the loop invariant hold. 
• (it < i'): Finally, consider an entry P[i t ,i'] = w a , where it < i' (consider an entry to 
the right of the diagonal). By the loop invariant, 

P[l,i'] = ^2 S S n { a , A i') w a , 
aeA' 

where A' C Ay is a subset of the set of entries above the diagonal in column i' . 
The row operation performed is sgn(Aj)(—l) St P[i t , i'] + P[l, i']. However, by Prop. 7, 
Aj U P[i t ,j] U P[it,i'] = A v . Therefore, w t G A v , and sgn(A 3 )(-l) St = sgn(a,Ai>). 
Moreover, we can be certain that a is not in the set A', since there is exactly one entry 
per column for each a = 1, . . . , n, and the only entries in A' are those that have been 
added there by previous row operations. Therefore, let A' U a = A" C A^ , and after 
the line 8 row operation, 

aeA" 

and condition (4) of the loop invariant holds. Moreover, when i' = i t + 1, A" = Ay, 
since every entry above the diagonal has been added to In this case, after i t 

is incremented to i t + 1 in line 9, condition (3) (the diagonal condition) of the loop 
invariant holds. Thus, when i t < i' , all relevant conditions of the loop invariant hold. 
Therefore, we have shown that if conditions (3) through (5) hold at the beginning of 
the j-th iteration, an arbitrary iteration of lines 6-9 will not alter those conditions. We 
have therefore shown that if conditions (1) through (6) hold at the j-th iteration, then 
conditions (1) through (6) also hold at the (j + l)-th iteration. 

Termination: 

The last row of Part Il [2™ _1 , :] consists of all zeros, except for the last n entries. In 
specific, Part„[2™ _1 , 2 n ~ 1 — n + 1 : 2™ -1 ]. Thus, during the iteration when column 
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j = 2" _1 — n + 1, the last row is added to the first row (with the correct signs since the 
loop invariant is maintained). Thus, during that iteration, P[l, :} = 0. Thus, when the 
algorithm terminates, the first row of the output matrix P consists entirely of zeros. This 
is the "useful information" that will help us prove our main result. □ 

Proof of Theorem 8 : The determinant of the partition matrix has degree 2 n_1 . We choose 
an arbitrary bitstring S G {0, l}™ -1 (which also defines an arbitrary partition of the set 
W), and substitute that factor into the partition matrix (line 2 of the algorithm Fac- 
torCheck. We have previously shown that the algorithm FactorCheck terminates 
with an output matrix P where the first row consists entirely of zeros. Therefore, the 
determinant of the output matrix P is zero, and the substitution performed in line 2 is 
a root of the determinant. Since the input factor defined by the bitstring S in the Fac- 
torCheck algorithm is an arbitary factor, and an arbitrary partition, we have shown 
that any factor S has the same result. Since every factor is a root, and there are 2"~ 1 
factors and only 2™ _1 roots of the determinant, the determinant of the partition matrix 
is the partition polynomial. □ 
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A. Appendix 

Here we display the FactorCheck algorithm running step-by-step input n = 4 and 
factor (— w\ — w 2 + w 3 + W4), described by bitstring S = 1100. Below follows P = Part4 
on line 1 of FactorCheck: 



P = 





123 


124 


134 


234 


1 


2 


3 


4 


1234 


W4 


W 3 


U)2 


Wl 














12 


W 3 


W4 








w 2 


Wl 








13 


W 2 










w 3 





Wl 





23 


w x 








W4 





w 3 


w 2 





14 





w 2 


w 3 





W4 








Wl 


24 





W\ 





w 3 





W4 





w 2 


34 








Wl 


W 2 








W4 


w 3 



















w 2 


w 3 


W4 



Below follows P on line 4, after the substition W4 = w\ + w 2 — w 3 , with column j = 1 
and row i = 2. During this iteration of the for loop, we perform row operations P[l, :] — 
P[2, :] + P[l, :], P[l, :] = -P[3, :] + P[l, :], and P[l, :] = -P[4, :] + P[l, :]. 





123 


124 


134 


234 


1 


2 


3 


4 


1234 


—^3 


w 3 


W'2 


Wl 














12 


W 3 


WI+W2—W3 








U>2 


Wl 








13 


UI2 


wi - 


(-102— 1«3 





U13 










23 


UIl 







-W2—W3 





">3 


t»2 





14 





W2 


WZ 













Wl 


24 





Wl 





UI3 


wi - 


(-102— 1«3 





W2 


34 








Wl 










UJ1+UJ2— WJ3 


W'i 





















W3 


wi+w 2 —w 3 



Below follows P on line 4 with column j '< = 2 and row i = 5. Here we perform the row 
operations P[l, :] = -P[5, :] + P[l, :] and P[l, :] = -P[6, :] + P[l, :]. 
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124 


134 


234 




1 










4 


1234 





Wl + W2 


— Wl + w 3 


-W2 + 


W3 


-W3 + 


W2 


Wl — W3 


— Wl — 


W2 





12 


W3 


W1+W2—VJ3 
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Wl 










13 


W2 





1^1+^2—^3 







W3 
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Wl 








W1+W2 - 


-W3 







W3 


W2 







14 





W2 


W3 







W1+W2 












Wl 


24 





Wl 





w 3 









W1+W2—VJ3 







W2 


34 








Wl 


W2 












W1+W2- 


-w 3 


W3 



















Wl 




W2 


W3 


wi- 


\-W2—u>3 



23 



Below follows P on line 4 with column j = 3 and row 
the for loop, we perform the single row operation P[l, :] 
both columns 3 and 4, row 1. 



= 7. During this iteration of 
P[7, :] :], which cancels 
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1234 
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— 11)2 


— tui — 
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— l«l — t«2 
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1«3 


U>2 







14 





U>2 
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24 
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W3 


i0i 


f 11i 2 —VJ3 







W2 


34 
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W2 








101+102 - 


-103 


W3 

















101 


UI 2 


W3 




101+102 —103 



Below follows P on line 4 with column j = 5 and row i = 8. During this iteration of the 
for loop, we perform the single row operation P[l, :} = P[7, :] + P[l, :], which cancels all 
the rest of the columns. 
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w 3 
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W2 


W3 


101+102—103 



Finally, below follows P upon termination. Notice that the first row of the matrix consists 
entirely of zeros, proving that this factor is a root of the determinant. 
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